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Variations in the broadband spectra of BL Lac objects: 
millimetre observations of an X-ray-selected sample 
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ABSTRACT 

Observations at millimetre wavelengths are presented for a representative sample of 22 
X-ray-selected BL Lac objects (XBLs). This sample comprises 19 High-energy cutoff 
BL Lacs objects (HBLs), 1 Low-energy cutoff BL Lac object (LBL) and 2 'inter- 
mediate' sources. Data for LBLs, which are mostly radio selected BL Lac objects 
(RBLs) are taken from the literature. It is shown that the radio-millimetre spectral 
indices of HBLs (q;5_23o = —0.29 ± 0.15) are slightly steeper than those of the LBLs 
(^5-230 — —0.19 ± 0.14). A correlation exists between a5_23o and 230 GHz luminos- 
ity. While this correlation could be an artefact of comparing two populations of BL 
Lac objects with intrinsically different radio properties it is also consistent with the 
predictions of existing unified schemes which relate BL Lac objects to Fanaroff-Riley 
class I radio galaxies. 

The HBLs have significantly flatter submillimetre-X-ray spectral indices 
(a23o-x = -0.68 ± 0.07) than the LBLs (a23o-x = -1.08 ± 0.06) although the 
two intermediate sources also have intermediate values of a23o-jf '^ —0.9. It is argued 
that this difference cannot be explained entirely by the viewing angle hypothesis and 
requires a difference in physical source parameters. The a2'io-x values for the HBLs 
are close to the canonical value found for large samples of radio sources and thus sug- 
gest that synchrotron radiation is the mechanism that produces the X-ray emission. 
As suggested by Padovani & Giommi, the inverse-Compton mechanism is likely to 
dominate in the LBLs requiring the synchrotron spectra of these sources to steepen 
or cutoff at lower frequencies than those of the HBLs. 

Key vifords: galaxies: active - BL Lacertae objects: general - radiation mechanisms: 
non-thermal - radio continuum: galaxies 



1 INTRODUCTION 

More BL Lac objects have now been discovered with X-ray 
surveys than with radio surveys. By definition, objects se- 
lected by the two techniques share many common properties 
such as compact core-jet structures, lack of or weak emission 
lines, rapid and large amplitude optical and radio variability, 
high and variable linear polarization and flat radio spectra. 
These general properties are universally attributed to the 
presence of a relativistic jet that is angled towards the line- 
of-sight (e.g. Blandford & Rees 1978). 

There are, however, several differences between XBLs 
and RBLs. Most strikingly, XBLs have significantly lower 
apparent radio luminosities than RBLs (e.g. Stocke et al. 
1990) whilst both classes have similar apparent X-ray lumi- 
nosities (Sambruna et al. 1994 and references therein). Ad- 
ditionally, XBLs typically have lower optical polarizations 
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and optical variability amplitudes (Schwartz et al. 1989; 
Jannuzi, Green & French 1993a; Jannuzi, Smith & Elston 
1993b), different parsec-scale radio properties (KoUgaard et 
al. 1996a), may exhibit negative evolution (Morris et al. 
1991; Perlman et al. 1996a; Bade et al. 1998) while RBLs 
exhibit positive evolution (Stickel et al. 1991) and have a 
larger galaxy fraction in their optical images (Wurtz 1994) 
and near-infrared colours (Gear 1993a). Furthermore, XBLs 
may exhibit steeper X-ray spectra than RBLs but this re- 
sult is only strictly true when comparing high-energy cutoff 
BL Lac objects (HBLs) which are mostly XBLs with low- 
energy cutoff BL Lac objects (LBLs) which are mostly RBLs 
(Padovani & Giommi 1996; Urry et al. 1996; see also Section 
2). 

Two competing models have been devised that attempt 
to explain at least some of these differences: (1) the beaming 
model in which XBLs are seen at larger angles to the line- 
of-sight than are RBLs (e.g. Ghisellini & Maraschi 1989; 
Padovani & Urry 1991; Ghisellini et al. 1993) and (2) a 
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model in which the cutoff frequency due to synchrotron ra- 
diative losses is at higher frequencies for XBLs (UV-X-ray) 
than for RBLs (optical-IR; Padovani & Giommi 1996 and 
references therein). 

For the latter model, the difference in radio power be- 
tween the two classes can be explained as a selection effect 
whereas the beaming model accounts for this difference in 
a natural way, although the jets have to accelerate between 
the X-ray emitting region and the radio emitting region in 
order to reproduce the observed behaviour. 

Recent studies suggest that the beaming model in com- 
bination with an intrinsic spread of physical jet properties 
may provide the best match to the available data (Sam- 
bruna, Maraschi & Urry 1996; Georganopoulos & Marscher 
1998; Fossati et al. 1998). 

This paper presents new millimetre wavelength obser- 
vations of a representative sample of 22 XBLs. A previous 
study that included data in this band (Gear 1993b) hinted 
that XBLs and RBLs exhibit similar radio-millimetre spec- 
tral indices but different millimetre-X-ray spectral indices, 
indicating that any differences between the two classes will 
be more apparent in the X-ray band. These results were 
based on a very small number of sources, however, and so 
the aim of this study is to provide a similar analysis for a 
larger number of objects. Throughout this paper we assume 
the Hubble constant. Ho — 50 Kms^^Mpc^^, Jlo ~ and 
spectral indices are written as 5*1^ ocv". 



2 SAMPLE 

BL Lac objects occupy a distinct region of the aox — Qro 
plane (Stocke et al. 1990). In fact, it was initially found 
that the RBLs and XBLs occupied two distinct regions 
but as more XBLs were discovered the separation became 
less pronounced (see below). The distribution is shown in 
Fig. 1 which includes RBLs from the '1 Jansky' sample 
(Stickel et al. 1991) and XBLs predominantly from the Ein- 
stein Medium Sensitivity Survey (EMSS; e.g. Perlman et 
al. 1996a), the Einstein Slew Survey (Perlman et al. 1996b) 
and the ROSAT all-sky survey (Kock et al. 1996; Nass et 
al. 1996; Bade et al 1998). Also included are sources from 
the Deep X-ray Radio Blazar Survey (DXRBS) which were 
found by correlating the ROSAT 'SNGkCAT database with 
several radio surveys (Perlman et al. 1998). Some BL Lac 
objects discovered in the ROSAT surveys tend to populate 
the aox — otro plane in the region between the '1 Jansky' 
RBLs and the Einstein XBLs and can be considered inter- 
mediate objects. Note also that many of the DXRBS sources 
have colours typical of the RBLs. 

We selected predominantly those sources with 4.85 GHz 
fluxes greater than 50 mjy but included fainter sources in 
an attempt to populate the aox — otro plane in a uniform 
manner. Those sources with detections at millimetre wave- 
lengths are shown as filled symbols in Fig. 1; millimetre XBL 
observations are either from this study (see Table 1) or from 
Gear (1993b) whilst the millimetre RBL observations are 
from Gear et al. (1994). 

Two of the sources, lES 0446-^449 and lES 0715-4259, 
have radio structures typical of radio galaxies (Perlman et 
al. 1996b) and are not considered as XBLs in the following 
analysis. Furthermore, it has been pointed out by Giommi 
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Figure 1. The distribution of BL Lac objects on the ciox — ciro 
plane. XBLs are represented by circles and '1 Jansky' RBLs are 
represented by triangles. Those sources highlighted with filled 
symbols have detections at 1.35 mm (~230 GHz). For complete- 
ness, new BL Lac objects from the DXRBS survey (see text) are 
shown as crosses; millimetre observations were not made for any 
of these sources. 



& Padovani (1994) and Padovani & Giommi (1995a) that 
the XBL/RBL distinction is not based on physical param- 
eters but only on the selection band and indeed several 
BL Lac objects are common to both samples. These au- 
thors proposed that BL Lac objects can be split into HBLs 
and LBLs based on the ratio of X-ray to radio flux density 
(fx/fr)- More recent results point to a continuous range 
of peak frequencies, making the HBL/LBL distinction it- 
self obsolete (e.g. Padovani 1999) but since this paper is 
largely concerned with a comparison of BL Lac objects se- 
lected from the two ends of the radio luminosity function 
the HBL/LBL approach is also considered here. The ra- 
tio log(/a;//r) (X-ray fluxes measured at 1 KeV and ra- 
dio fluxes at 5 GHz; both quantities in Jy) is calculated 
for both the '1 Jansky' sample and for the objects ob- 
served in Table 1. The LBL sample is defined as comprising 
those objects with iog{fx/fr) < —6.25, the HBL sample has 
^og{fx/fr) > —5.75 and an 'intermediate' sample, defined 
by —6.25 <log(/a;//r) < —5.75 contains only two sources, 
namely RXJ 12302-^2512 and EXO 1811.7-^3143. 

It turns out that, apart from the two intermediate 
sources, all of the XBLs considered in this paper are HBLs 
with the exception of 4C 47.08 which is a LBL and all of 
the RBLs are LBLs with the exception of Mrk501 which is 
a HBL. In Section 4 comparison is made between the LBL 
and HBL samples with comment on the relation of these to 
the intermediate sample but because of the large overlap, 
the same general results apply to the RBL/XBL samples. 
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Table 1. Summary of the observations 



source 


UT date 


2 mm 


1.35 mm 


0.85 mm o 


5-230 


a230-X 






(mjy) 


(mJy) 


(mJy) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




(7) 


RXJ 00079+4711 


19970921 


3cr < 17.3 
















19970624 




3cr < 16.3 




< 


-0.34 


> 


-0.84 




19980208 




3a < 7.5 




< 


-0.56 


> 


-0.78 




19980531 




3(T < 7.5 




< 


-0.56 


> 


-0.78 


lES 0033+595 


19970921 


34.3+8.1 


25.0+4.9 






-0.24 




-0.64 


4C 47.08 


19980207 




328.6+46.8 






-0.50 




-1.03 




19980531 


553.6+57.2 


395.9+43.5 






-0.45 




-1.04 


lES 0323+022 


19980309 
19980531 




13.6+2.9 
9.5+2.4 






-0.29 
-0.39 




-0.52 
-0.52 


lES 0414+009 


19980309 




24.2+3.6 






-0.28 




-0.63 


IBS 0446+449 


19970815 


3(T < 12.1 
















19970921 




3(T < 14.3 




<+0.18 


> 


-0.69 


lES 0502+675 


19980310 
19980311 




13.1+3.1 
18.2+3.3 






-0.24 
-0.15 




-0.63 
-0.65 


EXO 0556.4-3838 


19970921 


27.4+7.0 
















19980310 




3a < 12.1 




< 


-0.45 


> 


-0.48 


lES 0647+250 


19970815 
19970921 


36.3+8.5 


33.4+6.0 






-0.19 




-0.65 


lES 0715-259 


19970921 


3(T < 12.0 


3a < 6.5 




< 


-0.15 


> 


-0.65 


lES 0806+524 


19980310 




60.0+6.8 






-0.28 




-0.64 


IH 0829+089 


19970921 
19980310 


75.7+9.7 


3a < 8.4 
14.1+2.5 






-0.52 




-0.57 
-0.60 




19980424 


43.0+6.7 


22.9+3.8 






-0.40 




-0.63 


lES 1101-232 


19980505 


16.1+5.2 


14.9+4.2 






-0.39 




-0.57 


lES 1212+078 


19980505 
19980706 


3(7 < 27.4 


17.7+3.6 
16.7+3.3 






-0.44 
-0.45 




-0.76 
-0.76 


MS 1229.2+6430 


19980505 




12.6+4.0 






-0.31 




-0.64 




19980706 




3a < 11.8 




< 


-0.33 


> 


-0.64 


RXJ 12302+2512 


19980505 


225.4+27.3 


208.0+23.5 






-0.18 




-0.92 


EXO 1415.6+2557 


19980510 


3(T < 13.1 


3a < 12.3 




< 


-0.39 


> 


-0.58 


RXJ 16442+4546 


19970702 


17.7+4.4 


17.1+2.9 






-0.49 




-0.79 


4U 1722+11 


19970605 


118.8+23.7 
















19980207 


204.9+21.3 


188.2+19.3 


182.f 


±21.0 


+0.18 




-0.76 


IH 1730+500 


19980207 


57.2+9.0 


52.0+5.7 






-0.29 




-0.74 


EXO 1811.7+3143 


19980207 


46.4+8.3 


40.7+5.1 






-0.25 




-0.94 


lES 1959+650 


19970815 


144.0+19.3 


149.0+24.8 






-0.14 




-0.73 


MS 2143.4+0704 


19970815 
19980208 
19980531 
19980607 


14.3+4.7 
14.8+4.0 


3a < 21.8 
16.7+3.0 
11.8+2.5 




< 


-0.22 
-0.29 
-0.38 


> 


-0.79 
-0.77 
-0.74 


lES 2344+514 


19970624 
19970815 


78.9+11.0 
79.0+10.8 


54.1+10.4 






-0.36 




-0.74 



3 OBSERVATIONS 

The observations were made with the new submilhmetre 
camera, SCUBA (Gear et al. 1998, in preparation; Holland 
et al. 1998,1999), on the 15-m James Clerk MaxweU Tele- 
scope (JCMT), Mauna Kea, Hawaii. SCUBA, operating at 
75 mK, is optimised for mapping observations at submillime- 
tre wavelengths but also has single bolometers operating at 
2.0, 1.35 and 1.1 mm (the 1.1 mm detector was inoperable at 
the time of this work) corresponding approximately to 150, 
230 and 270 GHz respectively. The NEFDs at 2.0 and 1.35 
mm are typically measured as 120 and 60 mJyHz^^" re- 
spectively (including chopping) and these values are largely 
insensitive to the atmospheric opacity. Photometric obser- 
vations can also be made with the arrays, usually with the 



central pixel, and in this case the remaining bolometers can 
be used to subtract the sky emission which is correlated 
across the 2.3 arcmin field of view (Ivison et al. 1998; Jen- 
ness, Lightfoot & Holland 1998). Under average to poor 
conditions the NEFD at 0.85 mm (345 GHz) is ~ 90-200 
mJyHz^^/^. The standard SCUBA photometry mode was 
used for the observations. A small square map consisting of 
nine pixels and centred on the source is made with an inte- 
gration time of one second per pixel, the telescope then nods 
into the other beam where the procedure is repeated. The 
nine-point maps are then averaged to produce one photo- 
metric observation taking 18 seconds to complete. The pixel 
spacing was 3 arcsec at 2.0 and 1.35 mm and 2 arcsec at 
0.85 mm and the chop throw was 60 arcsec at ~7.8 Hz. 
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Figure 2. Broad-band spectra of the BL Lae objects. Millimetre data presented in this work are shown as filled symbols or upper limits. 



The observations were typically made under relatively 
poor observing conditions, i.e. high levels of precipitable wa- 
ter vapour (pwv) and atmospheric refraction (sky noise). 
Thus almost all data were taken at 2.0 and 1.35 mm to op- 
timise the sensitivity and reduce the effects of signal- loss due 
to sky noise and pointing errors (both giving variations of 
less than 3 arcsec whilst the beam sizes are in excess of 20 
arcsec at these wavelengths). The atmospheric opacity was 
monitored regularly with skydips at both 0.85 mm using 
SCUBA and at 1.3 mm with the nearby Caltech Submil- 
limeter Observatory (CSO) sky monitor. Standard relations 



were used to convert these opacities to those relevant to the 
adopted SCUBA filters. The instrumental gain was mea- 
sured on each night with observations of the planets Mars 
and/or Uranus. These gains are stable to within 10 per cent 
at 1.35 mm but may be elevation dependent at 2.0 mm with 
a total variation of order 20 per cent. We have assumed a 
calibration uncertainty of 10 per cent for all observations 
and these are added in quadrature with the uncertainty de- 
rived from the measured signal-to-noise ratio. Data reduc- 
tion was performed with the standard SCUBA reduction 
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Figure 2. - continued 



package and is analogous to that employed with previous 
instrumentation. 



4 RESULTS 

The observational results are summarized in Table 1 where 
column (1) has the source name, (2) the UT date of the ob- 
servation, (3), (4) and (5) the measured fluxes at 2.0, f.35 
and 0.85 mm, (6) the two-point spectral index, Q5-230 and 
(7) the two-point spectral index, a23o-x- See below for the 



definition of these quantities. Errors on the two-point spec- 
tral indices are not quoted since, because the data are non- 
simultaneous, it is impossible to quantify the effect of source 
variability. 

Fig. 2 shows the broad band spectra of those objects 
listed in Table 1. These spectra were constructed with non- 
simultaneous data taken from the literature but because of 
the broad frequency range considered, variability will have 
only a marginal effect on the overall spectral shapes. The 
two radio galaxies mentioned in Section 2 are included in 
this figure and it should be noted that the radio spectra in- 
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Figure 3. Histograms showing tiie 5 to 230 GHz spectral index 
(05-230) distribution for tiie HBLs (top) and LBLs (bottom). 



Figure 4. Histograms siiowing the 230 GHz to X-ray spectral 
index (0230- Jf) distribution for the HBLs (top) and LBLs (bot- 
tom). 



elude fluxes for both the compact and extended emission, the 
latter providing the dominant contribution. Neither source 
was detected at millimetre wavelengths. 

In general, the millimetre data fall close to a straight 
line connecting the 4.85 GHz point with the optical point 
and the millimetre flux is typically less than the radio flux 
as found for LBLs (Gear et al. 1994) and is thus likely to be 
optically thin unless the spectra have a complex structure. 

Columns (6) and (7) of Table 1 list the two-point spec- 
tral indices between 5 GHz (or more strictly 4.85 GHz) and 
230 GHz (Q5-230) and between 230 GHz and the X-ray band 
(q23o-x)- The radio and X-ray data, which are not contem- 
poraneous with the millimetre data, were taken from the lit- 
erature (Padovani & Giommi 1995b and reference therein). 
These spectral index distributions for the HBL sample are 
shown in the upper panels of Figs. 3 and 4 where two sources 
detected previously at 230 GHz (Gear 1993b) are also in- 
cluded, namely MS 1402.3-^0416 and MS 1534.2+0148 (both 
objects fulfill the HBL criterion defined in Section 2). 

The lower panels of these figures show the correspond- 
ing distributions for the LBL sample. The radio and mil- 
limetre data used to calculate the LBL spectral indices (and 
the HBL Mrk 501) were taken from Gear et al. (1994) and 
are quasi-simultaneous (taken within ~4 weeks) whilst the 
X-ray data (at 1 KeV) are from Urry et al. (1996) or, 
for 1514-241 and 2200+420 only, from Worrall & Wilkes 
(1990). Fig. 3 includes multiple epochs for individual sources 
since source variability can significantly affect the derived 
Q5-230 values even when the two fluxes are measured quasi- 
simultaneously (see Gear et al. 1994). In the case of non- 
detections, only the deepest upper limit is plotted. However, 
Fig. 4 only includes the flrst millimetre epoch of observation 
for each source taken from Gear et al. (1994) or Table 1. In 
this case, since variability will have little effect on 0230- x 
due to the large separation in frequency space, inclusion of 
multiple epochs would bias the statistics. 

The mean values of 05-230 are —0.29+0.15 for the HBLs 
and -0.19 + 0.14 for the LBLs. A Kolmogorov-Smirnov (KS) 
test shows that the Q5-230 distributions are different at the 
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Figure 5. Correlation of 05—230 with 230 GHz luminosity. Tri- 
angles are LBLs and circles are HBLs. The correlation coefficient, 
r, and its significance, P, arc indicated on the plot. 



99 per cent confidence level, or if the limits are assumed to be 
detections, 99.9 per cent confidence level, suggesting that the 
HBLs have slightly steeper spectra. Before concluding that 
this difference is an intrinsic property it is first necessary to 
investigate selection effects; for example, is Q5-230 correlated 
with either redshift or luminosity? 

Fig. 5 shows that 05-230 is indeed correlated with 
230 GHz luminosity in the sense that more luminous objects 
have more positive values of 05-230 (a Spearman rank-order 
correlation test gives r — 0.63 at the P => 99.99 per cent 
confidence level) whereas Fig. 6 shows no strong correlation 
between 05-230 and redshift (r = 0.32, P = 95.6 per cent). 
Subsequently, the correlation shown in Fig. 5 does not re- 
sult from a correlation of luminosity with both redshift (not 
shown but expected for flux limited samples) and spectral 
index. Possible origins of this correlation are discussed in 
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Figure 6. Correlation of 05-230 witli redshift. Triangles are 
LBLs and circles are HBLs. The correlation coefficient, r, and 
its significance, P, are indicated on the plot. 



Figure 8. Correlation of a230-x with redshift. Triangles are 
LBLs and circles are HBLs. The correlation cocfhcient, r, and its 
significance, P, are indicated on the plot. 
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Gear et al. 1993b). It is concluded that the millimetre to X- 
ray spectra of HBLs are significantly flatter than those of 
LBLs. 



5 DISCUSSION 

The preceding analysis showed that there are probable dil- 
ferences between the 05-230 distributions and distinct differ- 
ences between the a230-x distributions of LBLs and HBLs. 
This section discusses whether these results are consistent 
with a beaming model in which the LBLs are aligned closer 
to the line-of-sight than the HBLs or whether intrinsic dif- 
ferences between the classes are required by the data. 



5.1 The Q5- 



distribution 



Figure 7. Correlation of a230— X with X-ray luminosity. Trian- 
gles are LBLs and circles are HBLs. The correlation coefficient, 
r, and its significance, P, are indicated on the plot. 



the next section. It is thus found that the radio-millimetre 
spectral indices of LBLs are slightly steeper than those of 
HBLs. 

The Q230-X distributions (Fig. 4) are clearly different 
(>99.99 per cent confidence with a KS test). Mean values 
of Q230-X are -0.68 ± 0.07 for the HBLs and -1.08 ± 0.06 
for the LBLs. An important point to note is that the two 
intermediate sources that were discussed in Section 2 have 
intermediate values of a230-x (both fall in the —0.9 to —1.0 
bin; not shown). This point is discussed further in the next 
section. 

Figs 7 and 8 show the correlations of 0230 -x with X- 
ray luminosity and redshift. Fig. 7 shows no correlation (r = 
—0.05, P — 20.2 per cent) whereas Fig. 8 shows a possible 
anti-correlation of spectral index with redshift (r = —0.48, 
P = 98.85 per cent) which it most likely an artefact caused 
by the different redshift distribution of LBLs and HBLs (e.g. 



In Section 4 it was found that the HBLs have slightly steeper 
spectra between 5 and 230 GHz than the LBLs and that 
this result might follow from the correlation of Q5-230 with 
230 GHz luminosity. It is, however, possible that this corre- 
lation is spurious and results from comparing two groups of 
BL Lac objects with intrinsically different radio properties; 
note that Fig. 5 shows little overlap between the two clas- 
sifications. This section investigates possible origins for the 
correlation assuming it is real. 

One mechanism that can change both the luminosity 
and the spectrum is the relativistic Doppler effect, i.e. the 
beaming model. These models were first devised (Browne 
1983) to link LBLs with an unbeamed parent population of 
Fanaroff-Riley class I radio galaxies (FRI; Fanaroff & Riley 
1974). Subsequently, the HBLs, which display less extreme 
properties, were included in the same scheme but at larger 
viewing angles (less Doppler boosting) than the LBLs (Ghis- 
ellini & Maraschi 1989). Various studies suggest that the av- 
erage angles to the line-of-sight are ~ 10°, 20° and 60° for 
the LBLs, HBLs and FRIs respectively (Urry, Padovani & 
Stickel 1991; Ghisellini et al. 1993; KoUgaard et al. 1996b). 
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The details of the correlation shown in Fig. 5 will de- 
pend on the location of the synchrotron self-absorption turn- 
over frequency with respect to the frequencies used to calcu- 
late the spectral index. More specifically, a correlation will 
only be observed if the turn-over is Doppler shifted through 
the observed frequency range. 

The luminosity in the observer's frame, L^, is related to 
that in the source frame (denoted ') by Li, = 5''L'^ where 5 is 
the relativistic Doppler factor {5 — [Fp — (Fp — lY'^cosO]^^ 
where Fp is the pattern Lorentz factor and 6 is the view- 
ing angle). The exponent q = ?> — a for a spheroidal source 
whereas q = 2 — a for a continuous jet. The former is as- 
sumed below but note that the continuous jet case requires 
a higher Doppler factor to produce a given luminosity. Using 
the same notation, the observed frequency v — 5p' . 

The intrinsic synchrotron spectrum of a BL Lac object 
can be represented by 



L'^,-C'{.[)[^ 



exp 



a-i 



(1) 



where a is the optically thin spectral index, ^ is the spectral 
index in the optically thick regime, v[ is the frequency at 
which the optical depth equals unity and C'{vi) is a con- 
stant. It follows that 



L. 



''C'iv',) 



1 — cxp 



Sv[ 



Q-? 



(2) 



The intrinsic spectrum is constructed by assuming val- 
ues for a, ^, v'^ and C' {y[). Equation (2) can then be used to 
calculate the 230 and 5 GHz luminosities, and hence 05-230, 
in the observer's frame as a function of 5 (or S if a value for 
Fp is assumed) . Fig. 9 shows several of these simple models 
plotted over the observed correlation. The various parame- 
ters were chosen to give a good representation of the data 
and good agreement with the unified scheme, these being 
Q = -0.6, i = 0.4 and C'{v[) = 10^^ WHz-^; curves corre- 
sponding to values of i^i=3, 8, 12, 16 and 20 GHz are plotted 
from right to left. Various values of 5 are shown on the solid 
curve {v\ — 8 GHz). The corresponding values of 6 assume 
Fp = 5 which is the approximate value indicated by both the 
measured superluminal speeds (KoUgaard et al. 1996b, and 
references therein) and luminosity function studies (Urry & 
Padovani 1995). 

The model clearly provides a good representation of 
the observed correlation and is consistent with the unified 
scheme outlined above but are the assumed parameters rea- 
sonable? The optically thin spectral index, a — —0.6, is 
close to that measured for large samples of radio sources 
(e.g. Conway, Kellerman & Long 1963) and thus appears a 
valid assumption. As 8 approaches 90° the measured value 
of 05-230 will approach a as indicated in Fig. 9 and so, for 
example, reducing a will skew the envelope of models to the 
bottom right. Similarly increasing ^ will skew the envelope 
to the top left. 

The value of ^ is difficult to constrain from existing ob- 
servations. For a simple homogeneous source ^ — 2.5 but 
the composite spectra of BL Lac objects are well known to 
be flatter than this in the radio regime. The data for LBLs 
(Gear et al. 1994) give ^ in the range 0.1-0.9 but values 
towards both ends of this distribution give poor representa- 
tions of the observed correlation. 
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Figure 9. Simulated spectral index versus luminosity relation- 
ships plotted over the data shown in Fig. 5. The lines represent 
values of ui = 3,8, 12, 16,20 GHz (from right to left). Values of 
the relativistic Doppler factor, (5) and viewing angle {9 - see text) 
are shown on the solid curve for which ui = 8 GHz 



The value of C' {y[) fixes the various values of 5 and 
6 along the curves. Observations of FRI radio galaxies at 
230 GHz would show whether the assumed value is reason- 
able but no data exist in the literature at this time. However, 
estimates of 5 for BL Lac objects (mostly LBLs) exist in 
the literature from both Synchrotron Self-Compton (SSC) 
models (e.g. Ghisellini et al. 1993) and variability studies 
(Terasranta & Valtaoja 1994). The former study gives values 
of 5 = 0.01-14.3 (mean=3.9±4.2; 33 objects) with only four 
objects having 5 > 10 whilst the latter gives 5 — 0.2 — 4.3 
(mean=1.4± 1.2; 11 objects). These values are roughly con- 
sistent with those derived above and support the choice of 

It is thus concluded that the correlation of Q5-230 with 
230 GHz luminosity, if real, is consistent with that expected 
to arise if the LBLs and HBLs are Doppler boosted FRI 
radio galaxies. The parameters of the assumed model ap- 
pear reasonable but are perhaps limited to a narrow range 
of acceptable values. However, it is important to note that 
such a simple model cannot be used to infer that all of the 
differences between the spectra of BL Lac objects can be 
explained within the framework of the beaming model. For 
example, Sambruna et al. (1996) showed that whilst the ac- 
celerating jet model (Ghisellini & Maraschi 1989) at a range 
of viewing angles can reproduce successfully some of the 
spectral differences between LBLs and HBLs, it does not 
lead to the observed wide range of peak frequency or aro- 
These authors suggest that a difference in viewing angle cou- 
pled with an intrinsic difference in jet properties, such as 
magnetic field strength or jet size is required to explain the 
data. More recently, Fossati et al. (1997,1998; see also Sam- 
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bruna et al. 1996) have suggested that a hnk exists between 
spectral shape and bolometric luminosity and this idea has 
been extended to unify BL Lac objects in terms of viewing 
angle and jet electron kinetic luminosity (Georganopoulos 
& Marscher 1998). For the latter model, the difference in lu- 
minosity between HBLs and LBLs is only partly attributed 
to a difference in viewing angle. Furthermore, it can explain 
the differences in redshift distribution and cosmological evo- 
lution between the two classes but, as noted by the authors, 
does not include a treatment of the inverse- Compton mech- 
anism which may dominate the X-ray emission of LBLs (see 
next section). 



the synchrotron process but with a different electron en- 
ergy spectrum to the HBLs, although this situation appears 
somewhat contrived. A more realistic interpretation, taking 
into account the above discussion, is that the X-ray emission 
from the LBLs is inverse-Compton radiation. For the 'inter- 
mediate' sources it is natural to speculate that the X-ray 
emission comes from a combination of both mechanisms, re- 
quiring that the X-ray synchrotron emission is significantly 
steepened by radiative losses in these objects. 



5.2 The 0230 X distribution 

The spectral index distributions between 230 GHz and the 
X-ray band are significantly different and tend to cluster 
around values of ^ —0.70 for the HBLs and ~ —1.1 for 
the LBLs although, as noted in Section 4, the two 'interme- 
diate' BL Lac objects have intermediate values of a23o-x 
(~ —0.9). More millimetre observations of BL Lac objects 
that fall in the intermediate region of the aox — Qro diagram 
are needed to confirm this trend. The a23o-x distributions 
show a more pronounced difference between the two classes 
than the Orx distributions. For example, Sambruna et al. 
(1996) find a^^ = -0.57 ± 0.06 for the EMSS XBLs and 
ctrx = —0.86 ±0.08 for the '1 Jansky' RBLs with some over- 
lap between the two samples. This effect is likely to be due 
to the influence of synchrotron self-absorption on the radio 
fluxes whereas the 230 GHz emission is optically thin. 

The clear dichotomy in a23o-x certainly cannot be ex- 
plained within the framework of the simple beaming model 
outlined in the previous section but conversely cannot be 
used as evidence that beaming is unimportant in these 
sources. As discussed above the beaming model in combi- 
nation with one or more intrinsic differences between the 
two classes may provide the best flt to the data. 

The data presented here as well as in the literature 
now strongly suggest that different mechanisms are respon- 
sible for the X-ray emission in the two classes. Padovani 
& Giommi (1996) have analysed the ROSAT X-ray spec- 
tra of a large sample of BL Lac objects. They find that the 
LBLs are characterised by concave optical-X-ray continua 
and have flatter X-ray spectra than the HBLs. The HBLs 
have convex overall broadband continua and the X-ray spec- 
tral slopes are well correlated with aox and anti-correlated 
with the X-ray to radio flux and cutoff frequency. They con- 
clude that the X-ray emission in HBLs is produced by the 
synchrotron process whereas that from LBLs emanates from 
the SSC mechanism. In this scenario the HBLs must have 
high-energy cutoffs to the synchrotron emission at higher 
frequencies (UV-X-ray) than do the LBLs (IR-optical) oth- 
erwise synchrotron emission would dominate in both cases. 

The 230 GHz measurements presented here are on the 
high-frequency side of the self-absorption turn-over and pro- 
vide a measure of the optically thin spectral index. The mean 
value of — 0.68±0.07 for the HBLs is strikingly similar to the 
canonical value for large samples of radio sources (e.g. Con- 
way et al. 1963) providing strong evidence that the X-rays 
are produced by the synchrotron mechanism. 

For the LBLs there is no a priori reason to assume 
that the observed spectral indices are not also produced by 



6 SUMMARY 

Millimetre wavelength observations of 22 XBLs of which 19 
are HBLs, 1 is a LBL and 2 are 'intermediate' sources have 
been used along with similar published data for LBLs to 
search for differences between the radio-millimetre-X-ray 
spectral energy distributions of these objects. We flnd: 

(1) The Q5-230 two-point spectral indices of HBLs are 
slightly steeper than those of LBLs. Furthermore Q5-230 is 
correlated with 230 GHz luminosity. If this correlation is 
real, i.e. it is not a result of comparing two populations of 
BL Lac objects with intrinsically different radio properties, 
then it is consistent with existing unified schemes where BL 
Lac objects are FRI radio galaxies viewed along the jet axis, 
the HBLs being at larger viewing angles than the LBLs. 
The assumed parameters of this model, however, must take 
a relatively narrow range of values and it is unlikely that 
a difference in beaming parameters alone can unify LBLs, 
HBLs and FRIs. 

(2) The Q23o-jf two-point spectral indices of HBLs are 
significantly fiatter than those of LBLs and have values con- 
sistent with those expected for synchrotron radiation. The 
X-ray emission from the LBLs is more likely to be produced 
by the inverse-Compton mechanism. Interestingly, the two 
'intermediate' BL Lac objects also have intermediate values 
of a230-x suggesting that both processes contribute to the 
X-ray emission in these sources. 

The analysis presented in this paper thus points to a 
unified model for BL Lac objects that combines a spread of 
viewing angles with a spread of intrinsic source parameters. 
In this respect the model of Georganopoulos & Marscher 
(1998), adapted to incorporate inverse-Compton losses ap- 
pears promising. 
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